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Abstract 

A strengthened and generalized version of the standard Virtual Work Prin- 
ciple is shown to imply, in addition to bulk and boundary balances, a one- 
to-one correspondence between surface and edge hypertractions and hyper- 
strestress fields in second-grade continua. When edge hypertractions are con- 
stitutively taken null, the hyperstress is shown to take the form it has for a 
Navier-Stokes— a fluid, a relevant example of second-grade fluid-like material. 

1 Introduction 

The main conceptual point we want to make in this note is that stipulating a suit- 
able Principle of Virtual Powers to characterize mechanical equilibrium of continua 
of any grade bigger than one offers a key advantage: a Cauchy-type construction of 
the hyperstress fields accompanying the equilibrium hypertraction fields (a difficult 
task, that has been undertaken but not achieved so far) is no more needed, because 
hyperstresses can be explicitly computed in terms of hypertractions, and conversely. In 
this paper, we demonstrate this tenet in the case of second-gradient continua, by a 
simple argument. 

Theories of second-gradient continua have a long history. In the case of fluids, 
a possible dependence of pressure on the density gradient was first proposed by 
Korteweg Q9|] to model capillarity effects in 1901; for solids, two pioneering papers 
by Toupin HT5l 1 1611 on elastic materials with couple stresses appeared in the early 



1 



years 1960. A PVP approach to the formulation of the basic balance laws for these 
continua was taken by Germain J5j [6l in the early 1970s; the comprehensive 
article by Maugin 111 Oil appeared in 1980; a recent contribution, of special relevance 
to our present paper, is due to Gurtin and Fried O . 

It is well known that compatibility with the second law of thermodynamics for 
constitutive relations which are functions of the second deformation gradient de- 
mands that internal mechanical interactions have a nonstandard form. This has 
been discussed by many authors, mainly with an eye towards a better clarifica- 
tion of the role of the so-called hypertractions and hyperstresses; an important 
contribution was given by Dunn and Serrin [2]], who introduced the notion of 
the interstitial energy. An interesting feature of second-gradient materials is that, 
if bodies and subbodies having non everywhere smooth boundary are considered, 
then edge forces, that is, line distributions of hypertractions are to be expected (and, 
if a dependence on gradients higher than two is allowed, one has to deal also with 
vertex forces, as exemplified by Podio-Guidugli [12J). To our knowledge, a rigor- 
ous interaction theory accommodating such a nonstandard behavior remains to be 
constructed; interesting attempts in this direction have been carried out by Forte 
and Vianello O, Noll and Virga llllll . and Dell'Isola and Seppecher HTJ. 

Although here we do not deal with this difficult issue directly, in Section [3j the 
bulk of this paper, we do provide a full set of representation formulae not only, as 
is relatively easy, for tractions and hypertractions, both diffused and concentrated 
on edges, in terms of stresses and hyperstresses (see definitions (|25T )]. (1251) -,. and 
(1261 )) but also, conversely, for stresses and hyperstresses in terms of diffused and 
concentrated tractions and hypertractions (see d27l)- (l2"8"1) , and ([33])). Such represen- 
tation formulae generalize the corresponding formulae for simple (= first-gradient) 
materials, that we derive in our preparatory Section [2] Since we work in a non- 
variational setting, our results apply whatever the material response. The PVP we 
use includes edge tractions, both internal and external; without them, it would not 
be possible to arrive at the complete representation formula for the hyperstress in 
terms of hypertractions we construct in Subsection [33J 

Finally, in Section |4j we provide a new proof of the following not very well- 
known fact in the theory of second-gradient materials: if edge tractions are consti- 
tutively presumed null on whatever edge, then the hyperstress needs not be zero, 
although it takes a very special form whose information content is carried by a vec- 
tor field. We surmise that inability to develop edge interactions be characteristic 
of certain second-gradient fluids, an issue that we take up in a forthcoming paper 
II14II , continuing a line of thought proposed in 111 3 II . 

2 Simple Continua 

2.1 Power expenditures as constitutive requirements 

When a characterization of mechanical equilibrium is sought via a weak formula- 
tion of the virtual-work type, the primary object is a linear space "V of test(= virtual) 
velocity fields; the collection of tractions is introduced as the formal dual of "V, by 
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laying down a notion of external power expended in a virtual body motion; and 
the collection of stresses is introduced as the formal dual of the collection of test- 
velocity gradients, by laying down a notion of internal power. We regard specifica- 
tion of these two duality relations as the 'zeroth grade' of any constitutive theory. 

We classify a material body B as simple if the internal and external power expen- 
ditures have the following forms: 



T • gradv, 
tv, 



*T w (P)[v] := 

J 8P 

for all body parts(=subbodies) P and for all test-velocity fields v e "V . Here P is 
a bounded subset of the current observation space being regularly open (that is, 
coinciding with the interior of its closure) and having a part-wise regular boundary 
BP; we regard it as the region occupied by the typical body part at the time of 
our observation; time itself plays the role of a parameter. As to y, we accept the 
standard assumption that it includes all realizable velocities (that is, all velocity 
fields obtained by time differentiation of admissible deformation fields) and that it 
is closed under the operation of observer change, in the sense that, if v e y, then 

v+ = ^v) = q + Qv + Wx + 6f (1) 

for all translation velocities q, all rotations (= proper orthogonal tensors) Q, and 
all relative spins (= skew-symmetric tensors) W (for details, see the first subsection 
of the Appendix) . 

The internal field T, the Cauchy stress, is meant to measure the mechanical 
interaction of a material element of the subbody P with its immediate adjacencies; 
the external field t, the contact traction, is meant to account for the mechanical 
action exerted on P by its complement with respect to the material universe the 
body B belongs to. (We intentionally ignore all types of actions at a distance, 
because they are inessential to the purpose of our present discussion.) 

Both fields T and t are here introduced formally by way of Riesz duality with, 
respectively, the test velocity fields v and their gradient fields gradv. Both power 
expenditures T • gradv and t • v are required to be properly invariant under observer 
changes. Precisely, translational invariance of the external power expenditure over 
an arbitrary subbody implies that the contact traction field be balanced, i.e., that 



t = for all subbodies P. 

lap 

Moreover, the specific external power expenditure is rotationally invariant if and 
only if the contact traction is indifferent to observer changes, in the sense that 

0(t) = t + = Qt for all rotations Q. 

As to the specific internal power expenditure, which is quickly seen to be invari- 
ant under translational observer changes, its rotational invariance implies that, at 
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all points of B, the Cauchy stress field be symmetric-valued and indifferent to observer 
changes, in the sense that 

0{T) = T+ = QTQ r 
(for a proof of this result, see the Appendix) . 



2.2 Mutual consistency of stresses and tractions via a strength- 
ened Principle of Virtual Powers 

The mutual consistency of the stress and traction fields is the consequence of pos- 
tulating the following Principle of Virtual Power: 



W {i \P)[v]=W (e \P)[v], 



(2) 



for all body parts P and for all test-velocity fields v. 

Needless to say, the Principle is an invariant statement. The quantification on 
velocities is standard, that on body parts is not. Asking that ([2]) holds for all body 
parts is much stronger a requirement than demanding it to hold only for the whole 
bodyQ This additional strength connects the values taken by T and t at all points 
of B, and not only at its boundary points. 

With the use of a standard integration-by-parts lemma, we find that 



T-gradv: 



(-divT) • v + 



Tnv. 



CP 



Consequently, ([2]) can be written as follows: 



(-divT)-v-l- 



Tn- v 



dP 



t-V, 



dP 



for all body parts P and for all test velocity fields v; in particular, 



(-divT) -V + 



Tnv: 



SB 



t-v, 



(3) 



(4) 



dB 



for all test velocity fields v. Now, while both ([3]) and ([4]) imply the classic pointwise 
balance: 

- divT = at all points of B, 



it is ([3]) that implies that 

Tn = t at all points of B and for all unit vectors n, 



(5) 



and not only at the points of dB and for n the outward unit normal field, as implied 
by the weaker statement 



*We have been unable to assess who introduced this strenghtened qunatification in continuum me- 
chanics first, when and where. Needless to say, without it, it would not be possible to characterize 
equilibrium for a system of rigid bodies, nor the method of Euler cuts would make any sense. 
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A first direct consequence of §5$ is that 

t = t(x,n), 



namely, that - within the class of simple continua - the traction field at any point 
x e B must be thought of as depending on the orientation of the plane chosen 
through that point in order to detect the mutual mechanical contact interactions. 
Secondly it follows from ([5]) that the fields T and t carry essentially the same 
information. Indeed, while, given the Cauchy stress mapping x — > f (x), ([5]) yields 
that 



we also have that, given the traction mapping (x,n) <-* t(x,n), the Cauchy stress 
field can be constructed as follows: 



for n H (i = 1,2,3) any three mutually orthogonal directions. Thus, (O can be re- 
garded as a basic consistency condition for the pair of dynamic quantities (t, T) dual 
to the kinematic quantities (v, gradv), a condition that yields the representation 
formulae © and 0. 

All these results are more or less well known in the continuum mechanics com- 
munity. The main reason for us to recapitulate them is that they prompt the fol- 
lowing remarks, that are at the core of our present work. 

2.3 Thinking of complex continua 

Relations ([6]) and (0 are also arrived at when, as is customary only tractions on 
body parts are introduced, because stress is constructed d la Cauchy as a conse- 
quence of balance of tetrahedron-shaped parts. The Cauchy construction is the 
pillar on top of which the standard theory of diffuse (i.e., absolutely continuous 
with respect to the area measure) contact interactions stands. For complex (i.e., 
nonsimple) material bodies, a Cauchy-like construction has been attempted often, 
but not achieved so far, to our knowledge. Consequently, for such continua, al- 
though there is some sort of a general agreement about what notion of diffuse 
contact interactions applies, it is not clear what generalized stresses should ac- 
company them to achieve mechanical balance. Now, as we just showed for simple 
material bodies acted upon by diffuse contact interactions, a key advantage of a 
Principle of Virtual Power quantified over a large collection of body parts is that 
it dispenses us from going through Cauchy's tetrahedron construction to represent 
stress in terms of tractions. In the next section, this feature of a virtual-power ap- 
proach to formulate mechanical balance is demonstrated in the case of a popular 
type of complex continua: we show what consistency relations restrict the choices 
of tractions/hypertractions and stresses/hyperstresses, and we derive representa- 
tion formulae that generalize relations $6^ and (0, where edge hypertractions have 
a crucial role. 



i(x,n) = T(x)n, 



(6) 




(7) 
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3 Second-Gradient Continua 



3.1 Notation 

H, a capital letter from a sansserif font, is used to denote the third-order hyperstress 
tensor. We write H ij(t for the cartesian components of H with respect to a fixed 
orthonormal basis For any vector a, we write Ha for the second-order tensor 
with components (Ha) i; = Hij k a k (summation over repeated indexes understood); 
and, for any second-order tensor A, we write H [A] for the vector whose cartesian 
components are (H [A]); = Hij k Aj k ; in particular, H[a®b] = (Hb)a. 

Let y be the translation space associated to the euclidean space where 5f, a 
smooth surface, is embedded. Given any point of 5*, we denote by S P the orthog- 
onal projection of the vector space "V onto the tangent space to Sf at that point. 
The transpose of S P is % the inclusion map, which takes any vector in the tangent 
space into a vector of "V . Thus, for n a unit normal to the surface and I the identity 
mapping over "V, TP = I - n ® n@ 

3.2 Virtual power expenditures 

Among complex material bodies, second-gradient continua are those for which the 
internal and external power expenditures have the following forms: 

• (internal power expenditure) 



7T<0(P)[v] := 



r 



T • grad v + H ■ grad 2 v, (8) 



where grad 2 v := grad(gradv) is the second spatial gradient of the velocity 
field, and H is a third-order tensor field which satisfies: 

(Ha)b = (Hb)a, for all vectors a, b, (9) 

so as to have the same index symmetries as the field grad 2 v it is dual to; 

• [external power expenditure) 



^W(P)[v] 



(t-v + h- 3 n v) + 

dP 



_f e v, (10) 

dP 



where dP is the 'edgy' part of dP, if any (once again, and for the same 
reasons as before, any distance interaction entering the external power has 
been ignored). 

These definitions require some comments. Firstly, we here consider a part 
collection larger than usual: the boundary BP of a part may happen to be the 



2 We need not use here any of the well-known representations of S P and S I. The interested reader is 
referred to a paper by Gurtin and Murdoch Jsj . 



6 



union of a finite number of regular surfaces which join pairwise along edges, that 
is, smooth curves which begin and end at singular points called vertices; accord- 
ingly, the edgy part d P of d P is the union of a finite set of boundary curves where 
the surface-normal field has a jumpl^l Secondly, the introduction in ([8]) of the ad- 
ditional internal power expenditure associated with the hyperstress H is paraleled 
by the introduction in ( [To] ) of two types of external power expenditures associ- 
ated with hypertractions, both additional to the one associated with the traction 
field t, namely, the diffused hypertraction field h and the edge-force field f e . Thirdly, 
the consequences of requiring that both power expenditures be properly invariant 
under observer changes must be investigated. We treat the two last issues in the 
following subsections. 



3.2.1 Diffused and concentrated hypertractions 

Suppose that, for consistency with postulating an additional power expenditure in 
the bulk associated with the second velocity gradient, one associates with the first 
the following additional power expenditure at a part's boundary: 



H • grad v, 



?P 



with H some second-order tensor field work-conjugated to grad v. Then, on split- 
ting the velocity gradient into its tangential and normal parts: 

grad v = ( s grad v) s P + d n v <g> n, 



we have that 



H • grad v = 



SP 



(H 5 I)- s gradv + 



dP 



Hn-3„v. 



(11) 



3P 



Now, on setting Hn = h, the second integral on the right side of ( TTTb can be 
identified with the second surface integral in (flOl) : to motivate the presence of 
the line integral, we need a consequence of a surface-divergence identity, that we 
introduce right away. 

Let 5? be a regular surface, oriented by its normal n. At each point x of its 
boundary curve dy, a unit vector m can be chosen, orthogonal to both n and 
the tangent direction of dS? and pointing outward from the interior of 5^; such a 
vector m lies in the limiting tangent plane to y at x, and can be represented as 
m = s I/a, with /x the appropriate tangent vector. For any smooth tensor field A over 
Sf, the following integral identity holds true: 



(AT)- S gradv= _ 



sr 



s div(A s I)-v + 



5* 



Am v. 



(12) 



With a view to applying this identity to each regular portion of the boundary of 
an edgy body part, we draw attention to Figure[T] that depicts an edge cross-section 

3 A more precise and formal description of these geometrical notions is found in a paper by Noll & 
VirgaCIUl. 
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Figure 1: Four unit vectors describe an edge of a body part. 

of a part P. The limiting values of the outward unit normals to the surfaces y' and 
y" joining at the edge are n' and n", while the unit vectors m' and m" belong to 
the tangent planes to those surfaces, and point outward from their interior. The 
edge S is described by an ordered list of four unit vectors: 



S = (n^m'jn",!!!"). 



(13) 



Note that the vectors in the list are not independent, because they are coplanar 
and pairwise orthogonal; moreover, the ordered pairs (n',m') and (n",m") can 
be interchanged freely. With the stipulated conventions, ( TT3T) gives a local (first- 
order) description of the geometry of an edge, just as a unit normal provides a local 
(first-order) description of an oriented surface. 
With this notation, the identity (fl2l) yields: 



(H 5 I)- 5 gradv: 



JdP 



5 div(H s i)-v + 



dP 



[HmJ-v, 



(14) 



dP 



where [[HmJ denotes twice the edge average of Hm: 

[HmJ :=Hm' + Hm", 

a vector field over d P. Thus, while the first in fHl ) can be safely thought as ab- 
sorbed into the first surface integral in ( fTol l, the second motivates the introduction 
in (flOl) itself of the line integral over the edgy part of dP. 



3.2.2 Invariance of power expenditures under observer changes 

Translational and rotational invariances of the external power expenditure over an 
arbitrary subbody imply, respectively, the hypertraction balance 



t + 



f 6 = for all subbodies P. 



dp 



and the indifference properties 

t + = Qt, h + = Qh, (f e ) + = Qf e for all rotations Q. 
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It is shown in the Appendix that the rotational invariance of the specific internal 
power expenditure implies that, at all points of B, (i) the stress field T be symmetric- 
valued and indifferent to observer changes, just as it was the case for the Cauchy 
stress field in simple continua, and (ii) the hyperstress field H be indifferent to ob- 
server changes, in the sense that 

H+ = Q*H, with (Q*H)[a®b]-c=H[Q r a®Q r b]-Q r c. (15) 
3.3 The Principle of Virtual Powers and its consequences 

We postulate the following Principle of Virtual Powers: 

W ( - l \P)[v]=W ( - e \P)[v], (16) 

for all body parts P and for all test velocity fields v. 

We want to show that there are two types of consequences of the part-wise 
balance principle (fl6l ) : 

i. point-wise bulk and boundary balances, the same that would follow from a 
standard global Principle, stated for P = B only and quantified over the col- 
lection of test velocity fields; 

ii. mutual consistency relations for the stress fields (T, H) and the traction fields 
(t, h, f e ), such that either list determines uniquely the other, thus allowing 
the conclusion that each list conveys the same mechanical information. 

We begin by noticing that, via repeated integration by parts, we have that 



T-gradv+H-(grad 2 v): 



(T-divH)-gradv-l- 



Hn • gradv 



ep 



(-divT) -v + 



(Tn • v + ( H n)n • 3 n v + ( H n) • s grad v) , 



dP 



where we have set 



T := T - divH. 
Next, just as we deduced ( fT4l ) from ([T2l ). we find that 



(17) 



(Hn) - 5 gradv = - 



?p 



5 div((Hn)T)-v + 



dP 



CHn)m] - v. 



dP 



Thus, the PVP relation (fl6t can be given the following provisional form: 



(-divT) -v + 



+ 



J dP 

;Hn)m]-v: 



((Tn-Miv((Hn) 5 I)) • v+ (Hn)n • 3 n v) 



(t-v + h- 3 n v) + 



f e v, 



dP 



for all subbodies and all virtual velocity fields. 



(18) 
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By exploiting both quantifications with the use of smooth body parts only, we 
have from (PT8l ) that, at all points ofB, the following balances must hold: 

- divT = 0, (19) 

with T defined by fllTI ) in terms of T and H; moreover, for each unit vector n, 

Tn — s div((Hn) s I) = t, for all smooth surfaces oriented by n; (20) 

and 

(Hn)n = h. (21) 
Notice that the last relation can also be written as 

H [n <g> n] = h, for all unit vectors n. (22) 

The remain of (fT8l ) is 

I ([(Hn)m]-f e )-v=0, 

JdP 

for all edgy body parts and all test velocity fields. Consequently, we have that, again 
at all points ofB, 

[(Hn)m] =f e , for all edges S, (23) 

or rather, more explicitly, 

H [n' ® m' + n" ® m"] = f e , for all edges as in Figure Q] (24) 

Remark. Differences in notation apart, equations ( fT9l >, (f20l) -(f2Tb and (l23l coin- 
cide respectively with equations (4), (5) and (49) of Gurtin and Fried [4]; as they 
remark, when applied to the body's boundary, those equations are the same as 
equations (7.8), (7.9) and (7.10), derived variationally by Toupin in [15]. Our 
equation (1201) , however, is written in much more compact a form than the corre- 
sponding equations in |4]| and II15II . Interestingly, as noticed earlier by Forte and 
Vianello [3]], this opens the way to a suggestive analogy with 

Fn- s divM = t, 

the balance equation in shell theory that involves the force and moment tensors F and 
M, that is, the surface-stress descriptors of that theory; the following identifications 
suffice: 

F = T, M = (Hn) 5 I. 
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3.4 Bulk and boundary balances 

Relation JT9T ) is nothing but the 

• (bulk balance) 

-divT=0 inB. 

Relations ( 1201 ) and d2"TT) have, as we are going to see, a manifold use. 
Firstly, let (dB) r f es denote the portion of the body's regular boundary 

(3B) reg :=3B\3B 

where the external force and hyperforce fields t and h are assigned. Then, 
by a standard pill-box argument, we can deduce from ( 1201 ) the 

• [diffused traction & hypertraction boundary balances) 

Tn - s div((Hn) 5 I) = t , (Hn)n = h over (8B)^. 

Moreover, on (<3B) f , the portion of the edgy part of 8 B, if any, where external 
line forces f ^ are assigned, we deduce from ( [231) the 

• (concentrated hypertraction boundary balance) 

[(Hn)mfl=fe over( aB )f . 



3.5 Representation formulae 

Given the stress and hyperstress fields T and H over B U dB, relations ( fTTl ), ( 1201) , 
and (1211) . can be used to define the traction and hypertraction mappings: 

t^(x,n) := (T(x)-divH(x))n-Miv(H(x)n)% 
h(x,n) := (H(x)njn, 

for each unit vector n and all smooth surfaces 5? oriented by n; likewise, relation 
J24| ) can be used to define the line-traction mapping f e = f e (x, S): 

fe(x,#):=H(jc)[n / ®m / + ii"®m"], £ = (n',m';n",m"), (26) 

for all edges S, both of the body and of its parts. 

Our next goal is the deduction of a representation formula for H(x) in terms of 
the diffused and concentrated hypertractions acting, respectively, on three oriented 
coordinate planes through x and the relative coordinate edges; with such a formula 
at hand, it is easy to deduce a representation formula for T(x) that resembles (0): 

T(x) = 2(t(x,n I )+ 5 div(H(x)n, s I)) ®n„ (27) 
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n' = m" = e 2 




Figure 2: The coordinate edge & y 




Figure 3: The coordinate edges S 12 , S 23 , 8^ 



whence, for T(x): 

T(x) = T(x) + divH(x). (28) 

Remark. In addition to the complications that we are going to display along our 
path to determine H in terms of the surface and edge hypertraction mappings h 
and f e , the complicated dependence ( [271 ) of T on the diffused traction mapping 
t and on H itself gives a first idea of the difficulties intrinsic to any attempt to 
generalize Cauchy's tetrahedron construction we alluded at in Subsection 12.31 to 
second-gradient materials (let alone materials of higher grade) . 

At a given body point, let 5^ be the oriented plane parallel to e 2 ,e 3 whose nor- 
mal unit vector is e 1; and let the oriented planes 5^ 2 and y 3 be similarly defined. 
The coordinate edge S l2 is depicted in Figure [2j it is constructed as a special case 
of the edge in Figure [TJ by taking 5^ and Sf 2 for 5?" and Sf', respectively, and 
by choosing m' = e 2 , m" = e : (such prescriptions are needed to identify S 12 un- 
equivocally); the coordinate edges S 23 and S 13 are similarly defined (Figure PS]). 

With the use of (l22t and J24] l we define, respectively, the diffuse hypertractions 
h, on the coordinate planes 5^ and the concentrated forces f% on the coordinate 
edges S ik , namely, 

hj = !-![>! (Sex], 

h 2 = H[e 2 ®e 2 ], (29) 
h 3 = H[e 3 ®e 3 ], 
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and 

f e 12 = H[e 1 ®e 2 + e 2 8e 1 ]=2H[e 1 ®e 2 ], 

f e 13 = H[e 1 ®e 3 + e 3 ®e 1 ] =2H[e 1 ®e 3 ], 

f e 23 = H [e 2 ® e 3 + e 3 ® e 2 ] = 2 H [e 2 ® e 3 ], 

where we have accounted for the required symmetry of H, expressed by ([9]). We 
now show that the six vectors h; and f !v (each of which does depend on the Carte- 
sian basis chosen) allow for the construction of a basis-independent representation 
formula for the hyperstress tensor H . 
Indeed, we know that 

H = 2]Hy fc e i ®e J ®e fc ; 

ijk 

again taking the index symmetry of Hy fc into account, this sum can be rewritten as 
H = H w e, ® e ; « e, + 2 H >;fc e > ® e ; ® e * 

v-v- snsn , ^ (30) 

= 2j2j H w Ei ® e J ® e J + Zj 2j H ufc e i ® ( e ; ® e fc + e fc ® e ; ). 

j ' /<* ' 

Since 

Hijfc = e, -H[ey®e fc ] 

we deduce that 

1 

H ijk = ~ e. • H [ej ® e fc + e fc ® ej] . 

Thus, in view of ( 1291 ), 

1 1 

% = 2 e > ' H t e ; ® e J + e J ® e j] = 2 6 ' ' 2h J = 6 ' ' h J' 

which implies that 

hj = h(x, ej) = ^H i;j (x)e ; (no sum over j). (31) 

i 

Similarly for j < k, 

1 1 
H ijk = - e, ■ H [e, ® e k + e k ® e^] = -e ; • f 

which implies that 

f; fc = f e (x,^) = 2^]H IJ ,(x)e 1 , (forj<k). (32) 

i 

By substitution of ( f3TT ) and ( f32T ) in J30l ), we deduce the desired representation 
formula for H in terms of h and f e : 

H(x) = 2]h(x, e,) ® e ; ® e ; + - f e (x, <g} fe ) ® (e ; ® e fc + e k ® e,). (33) 
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We regard this formula as the main result of this paper. 



4 Second-gradient materials with no edge tractions 

We now ask the following question: which restrictions on H can be deduced from 
the assumption that no forces act on any edge, real or imagined? Said differently: 
is the hyperstress necessarily everywhere null in a second-gradient material body 
deemed to be unable to develop edge tractions? 

Let us then suppose that, at a fixed body point x, we have that f e (x, •) = 0, 
so that the forces on all coordinate edges are zero, whatever the basis. Under this 
assumption, we have from ( f33T ) that, with respect to any given basis, H can be 
writte as: 

H =^h J -®e J -®e i (34) 

;' 

(recall that vectors h ; do depend on the chosen basis). Now, consider a new or- 
thonormal basis e i; obtained from e ; by a rotation about e 3 of some angle 9: 

e x = cos 9 e 1 + sin 9 e 2 , e 2 = — sin 9 e x + cos 9 e 2 , e 3 = e 3 , 
and compute 

q 2 = 2H[e 1 ®e 2 ]. 

It turns out that 

i\ 2 = -2 sin0 cos <g> e-J +2 sin cos H[e 2 ®e 2 ] 

= h i =h 2 

+ cos 2 6 2H [ej ® e 2 ] - sin 2 9 2H [e 2 8 ej . 

=f e =f e 

But, because the edge tractions i\ 2 and t\ 2 must be zero by assumption, it follows 
from the last relation that, for all 9, 

= 2sin0cos0(h 2 -hj; 

this implies that hj = h 2 . Moreover, for a rotation about e 1 , the same argument 
yields that h 2 = h 3 . Then, we can set 

h := h : = h 2 = h 3 , 

and write ( |34l ) as follows: 

H = h ® y"| e ; ® e, = h ® I. 

i 

We conclude that a constitutive lack of edge forces implies that 

H=h®I, (35) 
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for some vector h; such h does not depend on the basis chosen for the argument, 
because 3h = H [I] . Conversely, if H admits the representation ( 1351) for some vec- 
tor h, then it readily follows from ( [241 and ( TT3T ) that all edge tractions are zero. 
Indeed, 

f e (*) = H[n'®m' + ii"®m"] 

= (l-(n / ®m' + n"®m"))h 
= (n'm' + n"m")h = 0, 
since (n',m') and (n",m") are pairs of orthogonal vectors. 

Remark. An alternative proof of the assertion leading to ( [351 was given by Forte 
and Vianello in Q3Q - Here it is. If there are no edge tractions, then 

w • H [n 8) m] = 

for any given vector w and for all pairs of orthonormal vectors n, m. Thus, the 
second-order symmetric tensor H T w, whose components are (H r w) jfc = w^,^), 
must have diagonal form with respect to any basis: 

wH = A(w)I, for some scalar-valued function A(w). 

However, since the left side of this relation is linear in w, there must be a vector h 
such that A(w) = h ■ w. But, 

wH = (h ■ w) I => H =h®I. 



Remark. No mathematically precise notion of fluidity and solidity has been sug- 
gested so far for complex continua. A current malpractice is to base the distinction 
on whether the first-gradient stress arising in response to deformation histories 
complies with the one or the other of the group-theoretic recipes proposed by Noll 
to sort simple fluids from simple solids. Recently, one of us [13] came up with a 
notion of macroscopic aggregation state of matter based on an invariance require- 
ment of the internal power expenditure under classes of changes of the reference 
placement different for fluid and solid simple materials. This notion is formally 
easy to generalize to complex materials, although the analysis involved is far from 
trivial. We develop these issues in full in a forthcoming paper [14]. Suffice it 
to mention here a class of Navier-Stokes—a fluids, that is, of second-gradient ma- 
terials obeying, among others to be introduced shortly, the following constitutive 
prescription for the power expenditure of hypertractions: 

H • grad 2 v = g ■ curl curl v . 

In view of the differential identity: 

curl curl v = grad(divv) — Av, 
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it is not difficult to see that 

g • curl curl v = H 1 ■ grad 2 v + H 2 ■ grad 2 v; 

with 

(Hi)yfc = -(5yg fc + 5 ifc ^) and H 2 = -g®I. 

Furthermore, if 

divv = 0, 

as is customarily stipulated in microfluidics, then 

g • curl curl v = — g • Av. 

Thus, (i) Hj can be interpreted as the reactive hyperstress associated to the second- 
gradient consequence of the standard incompressibility constraint; (ii) H 2 , the active 
hyperstress, induces no edge hypertractions; and, (iii) if one chooses the simplest 
constitutive prescription 

g=£Av, £> , 

then the regularizing term — ( A(Av) typical of the Navier-Stokes— a flow equation 
PI appears in the balance equation. 

5 Appendix 

We collect here some subsidiary material, with the purpose of improving the read- 
ability and self-containment of our paper. 

5.1 Observer changes 

By an observer change we mean the following transformation rule for position vec- 
tors of space points: 

x + = 0(x) = q + Qx, x := x - o, (36) 

where q is an arbitrary time-dependent translation and Q an arbitrary time-depen- 
dent rotation about point o. The related transformation rule for velocities is: 

v+ = q+Qv + Wx + , W:=QQ r = -W r , 

where all of the parameters q, Q, and W, can take arbitrary values at any chosen 
time (cf. CD). It follows from £36]) that 

grad+G) = grad(-)Q T . (37) 

Hence, as to the velocity gradient, one finds that 

(grad v) + = Q(grad v)Q r + W. (38) 
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5.2 Rotational invariance of the internal power expenditure 

The rotational invariance of the specific internal power expenditure is the require- 
ment that, at all points of B, 

T • grad v + H • grad 2 v = T+ • (grad v)+ + H + • (grad 2 v)+, (39) 

where 

T+ = 6(1), H + = 0(H), etc. 

In order to discuss the consequences of this requirement, it is convenient to de- 
fine the action of the rotation group on second- and third-order tensors: in cartesian 
components, 

[Q*A] y := Q ip Q jq A pq , [Q* H] yfc := Q ip Q ]q Q kr H pqr ; 

in absolute notation, 

Q*A=QAQ T , (Q*H)[a®b] -c= H[Q r a® Q T b] Q T c. 

Notice that, for any pair of rotations Q : and Q 2 , 

(Q1Q2) *A = Qj *(Q 2 * A), (Q1Q2)* H = Q X * (Q 2 * H); 

moreover, 

(Q*A)-B = A-(Q r *B), (Q*H)-K=H-(Q r *K), (40) 

for all pairs of second- and third-order tensors A, B, H, K and for all rotations Q. 
With this notation, the transformation rule ( 1381) becomes: 

(grad v) + = Q * (grad v) + W. 

Moreover, by an application of ( 1371 ), we readily deduce the transformation rule for 
the second gradient of velocity: 

(grad 2 v) + = Q * (grad 2 v). 

Thus, J39l can be rewritten as: 

T ■ gradv + H • grad 2 v = T + • (Q * grad v + W) + H + ■ (Q * grad 2 v), 

or rather, equivalently in view of ( 1401 ), 

T • gradv+ H ■ grad 2 v = (Q r * T + ) ■ gradv-l- T + ■ W + (Q r * H + ) ■ grad 2 v. 

Since all of the multipliers gradv, grad 2 v and W can be prescribed arbitrarily and 
independently from each other, we conclude that T must be symmetric and, more- 
over, that the following transformation rules must hold: 

T+ = Q*T, H + = Q*H 
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(cf. (flSt "). This conclusion fully generalizes the well known result for simple (first 
order) materials. 
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